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Abstract 

We give an overview over recent calculations of baryonic correlator functions 
with finite mass quarks in view on their applicability for QCD sum rules. The QCD 
sum rule method is then demonstrated within the Heavy Quark Effective Theory. 
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1 Mesonic and baryonic correlators 

Data sets which are produced to an huge amount especially by the so-called B factories 
like Barbar, Belle, Cleo and Hera B allow for the discovery of excited states of hadrons 
containing the h quark and other heavy quarks. In order to follow this development on the 
experimental side, theorists are asked to develop and use methods to analyze such excited 
states. There is one main obstacle on this way which consists of the scale differences 
that occur. While perturbative calculations including a heavy quark can effectively be 
done only in the high energy range, there is an extrapolation of these results down to the 
energy range needed where we expect excited states are to be located. 

One of the most powerful methods which was developed long ago are the QCD sum 
rules. In this talk we will only concentrate on the QCD sum rule approach developed 
by Shifman, Vainshtein and Zakharov, called SVZ approach []]]. This approach makes 
special assumptions about the spectral density of the correlator function related to the 
hadron and uses the Borel transform for the extrapolation. There are a lot of calculations 
for the correlator function of hadrons, especially for those containing a heavy quark. The 
calculations were performed within perturbative QCD as well as HQET and the massless 
limit. In the first part of this talk I will concentrate on calculations using perturbative 
QCD. In all cases the baryons are a kind of "stepchild" of the theorists, so we work hard 
to fill this gap. 

1.1 Baryonic correlators in QCD 

The mesonic correlator function for two vector currents has already been calculated ten 
years ago by Generalis 0] even in the case of two different and non-vanishing masses. 
Our aim is therefore to extend these calculations to the baryonic case. A first step has 
already been done by a recent publication ||. In this publication the spectral density is 
calculated for a three-quark current (two massless and one of finite mass) of the form 

Jb = [u lT Cd^ k e ljk . (1) 

The result presented in || is the one for the mass part IT m (— q 2 ) of the correlator 

II(-g 2 ) = ili q {-q 2 ) + mU m (-q 2 ). (2) 

This part is independently interesting because it can directly be compared with the result 
obtained within HQET, as we will see later. The momentum part U q (—q 2 ) is still under 
construction, we expect a publication in the next few months. As mentioned before, we 
can easily reconstruct the correlator function 

H m (-g 2 ) = r (3) 
if we know the spectral density. This is given by 

p-v-mpM- '=4°( 1+ H3H*' 1 } (4) 

where 

p (s) = 1 + 9z - 9z 3 - z 3 + 6z(l + z) In z (5) 
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and 

665 665 2 , /58 2 58 ,\, , 

Pl ( s ) = g + — z - —z 2 - g z 3 - ( — + 42z - 42z 2 - — z 3 ln(l - 2) + 
9 9 V 9 9 / 

+ (2 + l ^z - ^z 2 - ^z 3 ) \nz + A(l + 3z-3z 2 - ]-z 3 ) ln(l - z) \nz + 



9 J V3 3 ) 



+ 12z (2 + 3z + ^ 2 ) Q In 2 z - C(2)) + 4 (J + 12z + 3^ 2 - Li 2 (*) + 
+24,2(1 + z) (Uz{z) - C(3) - ^Li 2 (^) In . (6) 



This result is obtained by using basic integrals of the kind 

V(a, (3; g 2 /m 2 ) = J ^J )D {p2 + ^ _ p)W (?) 
which are a generalization of the standard object G(a,{3) of the massless calculation. 
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Figure 1: The ratio pi(s)/p (s) of the next-to-leading correction and the leading order 
term in dependence of the energy square s 



1.2 Comparison with the limits 

Here we only want to stress that we were able to compare this result with the two limits, 
i.e. the massless limit and the heavy quark or near-threshold limit. For the massless limit 
we obtain 

mp{s) = m m (p)s 2 |l + ^(21n (y) + f ) } (») 

where the relation between the pole (or invariant) mass parameter m and the MS mass 
m^g reads 




With this explicit finite mass representation we can even obtain terms like m 2 m(/z 2 /m 2 ) 
which are absent in the effective theory of massless quarks but result from perturbative 
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contributions of heavy quark condensates (^ r ^ r ). For the other limit, i.e. the near-threshold 
limit E — > with s = (m + E) 2 we obtain 

* - ^ + H£H ^ x + 4h &)) } + 0(£6) - (l0) 

The invariant function p m suffices to determine the complete leading HQET behaviour 
since one has <jip q + mp m — > m[i) + l)p for the leading term. In this region the appropriate 
device to compute the limit of the correlator is HQET (see e.g. @, H). Writing 

mp thT (m,E) = C(m/p,a s ) 2 p HQET (E,p) (11) 

we obtain the known result for p RClET (E, p) || with matching coefficient C(m/ p,a s ) @. 
In this case the matching procedure allows one to restore the near-threshold limit of the 
full correlator starting from the simpler effective theory near threshold || (see also 
Note that the higher order corrections in E/m to Eq. ( |T0"D can be easily obtained from 
the explicit result given in Eq. @. Indeed, the next-to-leading order correction in low 
energy expansion reads 

. thr/ _ 88E 6 ( a s (, ( p 2 \ 376 4vr 2 140, f m\\\ 
Ap (m, E) = -— |l + - (in ^ j + - + - + ^[-) j ) . (12) 

It is a much more difficult task to obtain this result starting from HQET. In Fig. p] 
we compare components of the baryonic spectral function in leading and next-to-leading 
order. Shown is the ratio pi(s)/p (s) where we put p, — m for simplicity. One can see 
that a simple interpolation between the two limits can give a rather good approximation 
for the next-to-leading order correction in the complete region of s. 




Figure 2: The calculated two- loop and three-loop topologies 



1.3 A specific feature of the calculation 

The two-loop and three-loop Feynman diagrams which had to be calculated for the deter- 
mination of the spectral function are shown in Fig. |2|. For our calculation we extensively 
used the fact that the diagrams can be composed by glueing together subdiagrams (see 
also Ref. [^D|). The so-called scalar spectacle diagram in Fig. ^ is calculated as the con- 
volution of two heavy-light spectral densities, 

p(s) = J A(s, si, s 2 )pi(s 1 )p2(s 2 )ds 1 ds2 (13) 

where the convolution function is given by the remaining line. 
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Figure 3: The scalar spectacle diagram 



2 QCD sum rules for the HQET 

The QCD sum rules for the correlator calculated above have not yet been constructed. 
Instead we will present the principles of the SVZ approach to QCD sum rules as applied to 
leading order HQET with results for heavy baryons taken from previous publications || 

0- 

2.1 Construction of QCD sum rules 

The QCD sum rules can be constructed by taking care of two possible expressions for the 
two-point correlator which in this case reduces to a scalar correlator function P(u>), 

n(u=p-v)=i J e vpx (T{J(x), J(0)})d 4 x = r'^^r , -Tr(rr)2 Tr(r^)P(u). (14) 

On the one hand side, the function P(co) satisfies the dispersion relation 

, f°° p(uj')duj' , 

p<u) = / Hy 1 + subtraction, (15) 

Jo uj' — to — iO 

where p(u) is the spectral density in HQET. On the phenomenological side the two-point 
correlator is represented by the spectral representation 

lip 12 lip 12 

P(uj) = T ^-^ + y 11 Xl + subtr. (16) 

A — CO — 20 UJx — — 2O 

where A = itlb — ttlq is the ground state energy of the baryon and Fb the residue. The 
main assumption of the SVZ approach is that the remaining sum can be approximated 
by the integral of the spectral density given by the dispersion relation and starting from 
some threshold energy Ec- The combination of the phenomenological and the theoretical 
identity for the correlator function then leads to 

2 ' ' + subtr. (17) 



A — us — iO Jo to' — uj — iO 
2.2 The Borel transformation 

This formula is not useful since the spectral density as calculated in the Euklidean domain 
is reliable only for negative values of u, while the integral is to be calculated mainly at 
u) — A. This region of integration can be reached by an extrapolation using higher and 
higher derivatives when uj goes to —00. This extrapolation is expressed by the Borel 
transformation (cf. e.g. Ref. 0) 

B^\f(uj))= lim ^ -T-/M, T = — fixed. (18) 
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The Borel transformation is by construction a derivative and therefore also cancels the 
(constant) subtraction terms. The Borel parameter T is an unphysical quantity in units 
of an energy, and the obtained values should be mostly independent on this parameter. 
This will be the main criterion in analyzing the sum rules. The Borel transformation 
leads to the final form of the QCD sum rule, 

±\F B \ 2 e- A / T = K(E c ,T), K(E c ,T) = J o ^e^du + Y,^. (19) 

At the end we can take the derivative with respect to the inverse Borel parameter and 
obtain a second sum rule, 




0.2 0.25 0.3 0.35 0.4 0.2 0.25 0.3 0.35 0.4 0.2 0.25 0.3 0.35 0.4 0.2 0.25 0.3 0.35 0.4 
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(a) (b) (c) (d) 

Figure 4: Sum rule analysis for the heavy baryon Aq. Shown are five curves for five 
different values of the threshold energy Ec spaced by 100 MeV around the central value 
E c = E c cst . E c grows from bottom to top. The figure parts are in detail (a) the lowest 
order sum rule results and (b) the 0(a s ) sum rule results for the bound state energy A, 
(c) the lowest order sum rule results and (d) the 0(a s ) sum rule results for the absolute 
value of the residue F A . 

2.3 QCD sum rule analysis 

The procedure of the sum rule analysis is as follows: First we calculate the theoretical 
expression for K(E C ,T) by calculating the spectral density within perturbation theory. 
This has been done for the next-to-leading order of the d = term as well as the d — 3 
term in the Operator Product Expansion. For brevity I will only mention the results of 
the analysis for the d — radiative corrections. For the analysis we have to select a "sum 
rule window", i.e. a range for the Borel parameter T in which the sum rule analysis is 
performed. The boundaries of this window is given by heuristic arguments. If the Borel 
parameter becomes too small, the nonperturbative contributions which are badly known 
blow up. If on the other hand the Borel parameter becomes too large, the fact that the 
Borel parameter appears as "temperature" in the sum rule causes the problem that higher 
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and higher excitations contribute to the ground state. The region of reliability and so the 
sum rule window is therefore roughly given by Aqcd < T < 2A. 

Now we use the second sum rule Q20D to determine the ground state energy parameter 
A. This is done by varying the continuum threshold parameter Eq in order to obtain a 
rather stable value for the quantities with respect to the unphysical parameter T within 
the sum rule window. The value obtained for the ground state energy can then be used in 
the first sum rule fll9D to determine the absolute value \F B \ of the residue. More detailed 
considerations are found in Ref. 0. Our results read 

E c (Aq) = 1.2±0.1GeV, 

A(A Q ) = 0.77 ± 0.05 GeV, (21) 
\F(A Q )\ = 0.022 ± 0.001 GeV 3 

while the 0{a s ) results are given by 

E c (Aq) = 1.1 ±0.1 GeV, 

A(A Q ) = 0.77 ±0.05 GeV, (22) 
\F(A Q )\ = 0.027 ±0.001 GeV 3 

where the errors are estimated by looking at the stability with respect to different values 
for Ec- Taking the experimental results for the masses of the baryons, namely m(A c ) = 



2284.9 ±0.6 MeV and m(A b ) = 5642 ±50 MeV gj, our central value A(A Q ) for the bound 



state energy suggests pole masses of m c = 1520 ± 100 MeV and m b = 4880 ± 100 MeV. 



3 Conclusion and Outlook 

The QCD sum rule analysis for heavy baryons requires the calculation of the spectral 
density related to the correlator function of baryonic currents. If this spectral density is 
calculated, the sum rule analysis can be used as powerful tool to extract phenomenological 
non-perturbative quantities like bound state energies and form factors. I have presented 
the sum rule analysis for an example within HQET. In order to perform a sum rule 
analysis for baryons containing finite mass quarks the following steps are still missing and 
will be done in the near future: 

• The calculation of the momentum part of the spectral density (nearly finished) 

• The matching procedure for the spectral density 

• The calculation of spectral densities for baryons with different quantum numbers 

Moreover it is planned to develop light-cone sum rules for the three-point function of 
heavy baryons in order to determine the Isgur-Wise function. 
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